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Abstract. The classical (scalar valued) theory of spherical functions 
(put forward by Cartan and others after him) allows one to unify un- 
der one roof a number of examples that were very well known before 
the theory was formulated. These examples include many special func- 
tions like Jacobi polynomials, Bessel functions, Laguerre polynomials, 
Hermite polynomials, Legendre functions, etc. 

All these functions had "proved themselves" as the work-horse in 
many areas of mathematical physics before the appearance of a unifying 
theory. Many of these functions have found interesting applications in 
signal processing in general as well as in very specific areas like medical 
imaging. 

The theory of matrix valued spherical functions gives a natural ex- 
tension of the well-known theory for the scalar valued case. 

The historical development in the matrix valued extension of this 
theory is entirely different. In this case the theory has gone ahead of 
the examples. 

The purpose of this note is to give some pointers to some examples 
and to "tease" some readers into this new aspect in the world of special 
functions. 

We close with a remark connecting the functions described here with 
the theory of matrix valued orthogonal polynomials. 

1. Introduction and statement of results 



The theory of matrix valued spherical functions, see [14] and [30], gives 



a natural extension of the well-known theory for the scalar valued case, see 



[ 23 1 . We start with a few remarks about the scalar valued case. 

The classical (scalar valued) theory of spherical functions (put forward by 
Cartan and others after him) allows one to unify under one roof a number of 
examples that were very well known before the theory was formulated. These 
examples include many special functions like Jacobi polynomials, Bessel 
functions, Laguerre polynomials, Hermite polynomials, Legendre functions, 
etc. 

All these functions had "proved themselves" as the work-horse in many 
areas of mathematical physics before the appearance of a unifying theory. 
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Many of these functions have found interesting applications in signal pro- 
cessing in general as well as in very specific areas like medical imaging. It 
suffices to recall, for instance, that Cormack's approach || — for which he 
got the 1979 Nobel Prize in Medicine, along with G. Hounsfield — was based 
on classical orthogonal polynomials and that the work of Hammaker and 



Solmon [20 1 as well as that of Logan and Shepp |2(| is based on the use of 
Chebychev polynomials. 

The crucial property here is the fact that these functions satisfy the inte- 
gral equation that characterizes spherical functions of a homogeneous space. 
For a review on some of these topics the reader can either look at some 
of the specialized books on the subject such as [23] or start from a more 
introductory approach as that given in either |8| and 

This integral equation is actually satisfied by all Gegenbauer polynomi- 
als and not only those corresponding to symmetric spaces. This point is 
fully exploited in Q where this property is put to use to show that different 
weight functions can be used in carrying out the usual tomographic oper- 
ations of projection and backprojection. This works well for parallel beam 
tomography but has never been made to work for fan beam tomography be- 
cause of a lack of an underlying group theoretical formulation in this case. 
For a number of issues in this area, including a number of open problems, 
see @. 

For a variety of other applications of spherical functions one can look at 

i , m, & 

We now come to the main issue in this note. 

The situation with the matrix valued extension of this theory is entirely 
different. In this case the theory has gone ahead of the examples and, in 
fact, to the best of our knowledge, the first examples involving non-scalar 
matrices have been given recently in [11] , [|l^] and [13]. For scalar valued 



instances of non-trivial type, see [22 



The issue of how useful these functions may turn out to be as a tool in 
areas like geometry, mathematical physics, or signal processing in the broad 
sense is still open. From a historical perspective one could argue, rather 
tautologically, that the usefulness of the classical spherical functions rests 
on the many interesting properties they all share. With that goal in mind, it 
is natural to try to give a glimpse at these new objects and to illustrate some 
of their properties. The rather mixed character of the audience attending 
these lectures gives us an extra incentive to make this material accessible to 
people that might normally not look in the specialized literature. 

The purpose of this contribution is thus to present very briefly the essen- 
tials of the theory and to describe one example in some detail. This is not 
the appropriate place for a complete description, and we refer the interested 



reader to the papers [11] , [T^] and [13]. 

We hope to pique the curiosity of some readers by exploring the extent 
to which the property of "positive linearization of products" holds in the 
case of the spherical functions associated to P2(C). This result has been 



MATRIX VALUED SPHERICAL FUNCTIONS 



3 



important in the scalar case, including its use in the proof of the Bieberbach 
conjecture, see [^]. The property in question is illustrated well by considering 
the case of Legendre polynomials: the product of any two such is expressed 
as a linear combination involving other Legendre polynomials with degrees 
ranging from the absolute value of the difference to the sum of the degrees 
of the two factors involved. Moreover, the coefficients in this expansion are 
positive. 

We should stress that the intriguing property described here is one enjoyed 
by a matrix valued function put together from different spherical functions 
of a given type. In the classical scalar valued case these two notions agree 
and the warning is not needed. This combination of spherical functions has 



already been seen, see [11] , [12] and [13] to enjoy a natural form of the 



bispectral property. For an introduction to this expanding subject we could 
consult, for instance, [16], [17]. The roots of this problem are too long to 
trace in this short paper, but the reader may want to take a look at [15]. 
For off-shoots that have yet to be explored further one can also see |i8}] and 



[19|. The short version of the story is that some remarkably useful algebraic 
properties that have surfaced first in signal processing and which one would 
like to extend and better understand have a long series of connections with 
other parts of mathematics. For a collection of problems arising in this area 
see 

The issue of linearization of products, without insisting on any positivity 
results, plays (in the scalar valued case) an important role in fairly success- 
ful applications of mathematics. For example, the issue of expressing the 
product of spherical harmonics of different degrees as a sum of spherical har- 
monics plays a substantial role in both theoretical and practical algorithms 
for the harmonic analysis of functions on the sphere. For some developments 



in this area see p] as well as [25 



In the context of Quantum Mechanics this discussion is the backbone of 
the addition rule for angular momenta as can be seen in any textbook on 
the subject. 

In the last section we make a brief remark connecting the functions de- 
scribed here with the theory of matrix valued orthogonal polynomials, as 
developed for instance in || and 0. 



2. Matrix valued spherical functions 

Let G be a locally compact unimodular group and let K be a compact 
subgroup of G. Let K denote the set of all equivalence classes of complex 
finite dimensional irreducible representations of K; for each 5 £ K, let £5 
denote the character of 5, d(8) the degree of 5, i.e. the dimension of any 
representation in the class 5, and xs = d(5)£s- 
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Given a homogeneous space G/K a zonal spherical function (|23[]) ip on 
G is a continuous complex valued function which satisfies (/9(e) = 1 and 

(1) (p(x)<p(y) = (p(xky)dk x,y£G. 

Jk 

A fruitful generalization of the above concept is given in the following 
definition. 



Definition 2.1 (|p0fl,p4|). A spherical function on G of type 6 G K is a 
continuous function on G with values in End(F) such that 

i) <E>(e) = / (i= identity transformation). 

ii) ®(x)®(y) = j K xs{k~ l )^{xky)dk, for all x,y G G. 

The connection with differential equations of the group G comes about 
from the property below. 

Let D(G) K denote the algebra of all left invariant differential operators 
on G which are also invariant under all right translation by elements in 
K. If (V, 7r) is a finite dimensional irreducible representation of K in the 
equivalence class 5 G K, a spherical function on G of type 5 is characterized 
by 

i) $ : G — > End(V) is analytic. 

ii) ®(kigk 2 ) = ir(ki)®(g)ir(k 2 ), for all k\, k 2 G K, g G G, and $(e) = /. 

iii) [D$](g) = $(g)[D$](e), for all D G D(G) K , g G G. 

We will be interested in the specific example given by the complex projective 
plane. This can be realized as the homogeneous space G/K, where G = 
SU(3) and K = S(U(2) x U(l)). In this case iii) above can be replaced 
by: [A 2 &](g) = A 2 <I ) (s'), [A 3 $](#) = \s^(g) for all g G G and for some 
A2, A3 G C. Here A2 and A3 are two algebraically independent generators of 
the polynomial algebra D(G) of all differential operators on G which are 
invariant under left and right multiplication by elements in G. 

The set K can be identified with the set Z x Z>o- If k = (q °), with 
A G U(2) and a = (det^)" 1 , then 

7r(fc) =KnM) = (detA) n A 1 , 

where A 1 denotes the /-symmetric power of A, defines an irreducible repre- 
sentation of K in the class (n,l) eZx Z>o- 

For simplicity we restrict ourselves in this brief presentation to the case 
n > 0. The paper [11] deals with the general case. The representation 7T n j 



of U(2) extends to a unique holomorphic multiplicative map of M(2, C) into 
End(V^), which we shall still denote by 7r n ;. For any g G M(3,C), we shall 
denote by A(g) the left upper 2x2 block of g, i.e. 



A{g) 



gn 912 

,921 922, 

For any ir = 7T( n n with n > let & n : G — > End(V^-) be defined by 

MS) = ®nM = KnMi.9))- 
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It happens that ^ is a spherical function of type (n, I), one that will play 
a very important role in the construction of all the remaining spherical 
functions of the same type. 
Consider the open set 

A = { g G G : det A(g) ^ } . 

The group G = SU(3) acts in a natural way in the complex projective plane 
P 2 (C). This action is transitive and K is the isotropy subgroup of the point 
(0,0,1) G P 2 (C). Therefore P 2 (C) = G/PC. We shall identify the complex 
plane C 2 with the affine plane { (x, y, 1) G P 2 (C) : (x, y) G C 2 }. 

The canonical projection p : G — ► P2(C) maps the open dense subset 
A onto the affine plane C 2 . Observe that „4 is stable by left and right 
multiplication by elements in K. 

To determine all spherical functions : G — ► End(V n -) of type n = ir nt i, 
we use the function $ n introduced above in the following way: in the open 
set A we define a function H by 

H(g) = $(g) ^{g)~\ 

where <3? is suppose to be a spherical function of type ir. Then H satisfies 

i) H(e) = I. 

ii) H(gk) = H(g), for all g G A, k G K. 

iii) P(%) = K{k)H{g)K{k~ l ), for all G A k G 

Property ii) says that H may be considered as a function on C 2 . 

The fact that <I> is an eigenfunction of A 2 and A3 makes H into an eigen- 
function of certain differential operators D and i? on C 2 . 

We are interested in considering the differential operators D and E applied 
to a function H G C°°(C 2 ) ® End(K) such that H(kp) = ^(fc)F(p)^(A;)- 1 , 
for all k G K and p in the affine complex plane C 2 . This property of H 
allows us to find ordinary differential operators D and E defined on the 
interval (0, 00) such that 

(DH)(r,0) = (DH)(r), (M)(r,0) = (M)(r), 

where H(r) = H(r,0). 

Introduce the variable t = (1 + r 2 )~ 1 which converts the operators D and 
E into new operators D and E. 

The functions H turn out to be diagonalizable. Thus, in an appropriate 
basis of V w , we can write H(r) = H(t) = (ho(t), ■ ■ ■ , hi(t)). 

We find it very convenient to introduce two integer parameters w, k sub- 
ject to the following three inequalities: < w, < k < I, which give 
a very convenient parametrization of the irreducible spherical functions of 
type (n,l). In fact, for each pair (l,n), there are a total of / + 1 families 
of matrix valued functions of t and w. In this instance these matrices are 
diagonal and one can put these diagonals together into a full matrix valued 
function as we will do in the next two sections. It appears that this function, 



F. A. GRUNBAUM, I. PACHARONI, AND J. TIRAO 



which concides with the usual spherical function in the scalar case, enjoys 
some interesting properties. 



The reader can consult [11] to find a fairly detailed description of the 
entries that make up the matrices mentioned up to now. A flavor of the 
results is given by the following statement. 

For a given I > 0, the spherical functions corresponding to the pair (l,n) 
have components that are expressed in terms of generalized hypergeometric 



functions of the form p _|_2-F p +i, namely 



3=0 



3. The bispectral property 

For given non-negative integers n, I and w consider the matrix whose rows 
are given by the vectors H (t) corresponding to the values k = 0, 1, 2, . . . , I 
discussed above. Denote the corresponding matrix by 

$(w,t). 

As a function of t, &(w,t) satisfies two differential equations 

D$(w, tf = $(w, t)*A, E$>{w, tf = $(w, t)*M . 

Here A and M are diagonal matrices with A(i, i) = —w(w + n + i + l + \) — 
(i - l)(n + i), M(i, i) = A(i, i)(n - I + 3i - 3) - 3(i - - i + 2)(n + i), 
1 < i < I + 1; D and E are the differential operators introduced earlier. 
Moreover we have 

Theorem 3.1. There exist matrices A w , B w , C w , independent oft, such 
that 

Au,$(w -l,t) + B w ®(w, t) + C w ®(w + 1, t) = t$(w, t) . 

The matrices A w and C w consist of two diagonals each and B w is tridi- 
agonal. Assume, for convenience, that these vectors are normalized in such 
a way that for t = 1 the matrix &(w, 1) consists of all ones. 

For details on these matrices as well as for a full proof of this statement, 



which was conjectured in [11], the reader can consult 



4. Linearization of products 

The property in question states that the product of members of certain 
families of (scalar valued) orthogonal polynomials is given by an expansion 
of the form 

j+i 

k=\j-i\ 

and that the coefficients in the expansion are all non-negative. 
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For a nice and detailed account of the situation in the scalar case, see for 
instance fl|, [p7| . Very important contributions on these and related matters 
are |J and |24j| . 

It is important to note that the property in question is not true for all 
families of orthogonal polynomials, in fact it is not even true for all Jacobi 
polynomials P^*'^\ normalized by the condition Pi, a '^(l) positive. For our 
purpose it is important to recall that non-negativity is satisfied if a > (3 and 
a + p > 1. 

The case I = 0, n > 1. 

From |TI[ we know that when I = and n > the appropriate eigenfunc- 
tions (without the standard normalization) are given by 

$( W ,t) = 2 F 1 (~ w >Xi n+2 > t )- 

This means that with the usual convention that the Jacobi polynomials are 
positive for t = 1 we are dealing with the family 

P^\t). 

If n = or n = 1 the family Pw' n ^ meets the sufficient conditions for non- 
negativity given above. For n = the coefficients a>k are all strictly positive; 
in the case n = 1 the coefficients a\i-j\+2ki are strictly positive while the 
coefficients a\ i _j\ + ! e , k odd, are zero, as the example below illustrates. 

We now turn our attention to the case n > 1. 

Conjecture 4.1. For n an integer larger than one, the coefficients in the 
expansion for the product PiPj above alternate in sign. 

This conjecture is backed up by extensive experiments, one of which is 
shown below. It deals with the case of w (i.e. i and j ) equal to 3 and 4. 
Prof. Richard Askey, supplied a proof of this conjecture. This gives us a 
new chance to thank him for many years of encouragement and help. 

The product of the (scalar valued, and properly normalized) functions 
<£(3,t) and $(4, £) is given by the expansion 

$(3, t)$(4, t) = oi*(l, t) + a 2 $(2, t) + a 3 $(3, t) + a 4 $(4, t) 
+ a 5 $(5, t) + a 6 $(6, t) + a 7 $(7, t) 

with coefficients given by the expressions 
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ai 
a 2 



(n+2)(ra+3)(n+4) 
(n+8)(n+9)(n+10) 

6(n- 1) (n+3)(n+4) (n+6) 2 
(n+7)(n+8)(n+9)(n+10)(n+ll) 

3(ra+4)(n+5)(7ra 3 +52n 2 +67n+162) 



0,3 (n+7)(Vi+9)(n+10)(n+ll)(n+12) 

4(n-l)(n+6)(llra 3 +123Tt 2 +436ra+648) 
(n+8)(n+9)(n+ll)(n+12)(n+13) 

3(ra+5)(ra+6) (n+7) (19n 3 +155n 2 +162n+504) 
(n+8)(n+9)(n+10)(n+ll)(n+13)(n+14) 

42(n-l)(n+5)(n+6) 2 (n+7)(n+8) _ 
(n+9)(n+10)(n+ll)(n+12)(n+13)(n+15) 

14(n+5) (n+6) 2 (n+7) 2 (n+8) 



04 

a 5 
a 6 



0,7 (n+10)(n+ll)(n+12)(n+13)(n+14)(n+15) 

This shows that even in the scalar valued case, as soon as we are dealing 
with non-classical spherical functions we encounter an interesting sign al- 
ternating property that is quite different from the more familiar case. Here 
and below we see that things become different once n is an integer larger 
than one. 

Now we explore the picture in the case of general I . 
The case I > 0, n > 1. 

Conjecture 4.2. If i < j then the product of &(i,t) and ®(j,t) allows for 
a (unique) expansion of the form 

j+i+l 

$(i,t)$(j>*) = Yl A kHk,t). 

fc=min{j— i— 1,0} 

Here the coefficients are matrices and the matrix valued function <&(w,t) 
is the one introduced in section 3. This conjecture holds for all nonnegative 
n and is well known for I = and n = 0. 

One should remark that in the case of I = we obtain the usual range 
in the expansion coefficients ranging from j — i to j + i as in the case of 
addition of angular momenta. For larger values of I we see that extra terms 
appear at each end of the expansion. 

Conjecture 4.3. If i < j then the coefficients in the expansion 

j+i+l 

fc=min{j— i— 1,0} 

with k in the range j — i,j + i have what we propose to call "the hook 
alternating property. " 

We will explain this conjecture by displaying one example. First notice 
that we exclude those coefficients that are not in the traditional or usual 
range discussed above. 
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At this point it may be appropriate in the name of truth in advertisement 
to admit that we have no concrete evidence of the significance of the property 
alluded to above and displayed towards the end of the paper. We trust that 
the reader will find the property cute and intriguing. It would be very 
dissapointing if nobody were to find some use for it. 

The results illustrated below have been checked for many values of I > 0, 
but are displayed here for I = 1 only. 



Recall that from [11] the rows that make up the matrix valued function 



H (t, w) are given as follows: the first row is obtained from the column vector 
H(t) = 



1 A \ r / —IDs w+n+3, A— n 
L ~TiTl) 3t2 \ n+2, A-n-1 ' 



I —w, w+n+3 , 
VI ( n+1 >* 

with 

A = —w(w + n + 3) 

and the second row comes from the column vector 

-w, ui+n+4 

/ VI I 

H(t) 



2* 1 { n+2 



I i -s \ T7i I —w—1, w+n+3, A 

-(n + l)^ n+1)A _! ;t 
with 

A = —w(w + n + 4) — n — 2. 
The product of the matrices <£(2,i) and <£(6,i) is given by the expansion 

$(2, t)*(6, t) = A 3 $(3, t) + A 4 $(4, t) + A 5 $(5, t) + A 6 $(6, t) 
+ A 7 $(7, t) + A 8 $(8, t) + A 9 $(9, t) 
where the coefficient matrices are displayed below. The matrix 
[0 \ 

^3 = n 16(n+4)(n+5)(n+6) 2 (n+7) 2 • 

\ U (n+11) (n+12) (n+13) (n+14) (n+15) (n+16) / 



The matrix 



A 4 



L21 L22 



with 
L 



15(n+5) 2 (n+6)(n+8) 
11 ~~ 2(n+12) (n+13) (n+14) (n+15) ' 



j _ 5(n+5)(n+6)(4 n 2 +55n+216) j (n+5)(n+6)(n+7)(8 n 2 +153n+724) 

12 ~~ 6(n+13) (n+14) (n+15) (n+16) ' 21 ~~ 2(n+12) (n+13) (n+14) (n+15) (n+16) ' 

j __ 5(n+6)(n+7)(248 n 4 +4665 ra 3 +27202 n 2 +45137 n-23252) 
^ 22 ~~ 12(n+ll)(n+13)(n+14)(n+15)(n+16)(n+17) - 

The matrix 

_ M n M 12 
M 21 M 22 
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with 

M (n+5)(n+6)(185 n 3 +3284 n 2 +15732 n+10368) 
iWn ~~ 6(n+7)(n+12)(n+14)(n+15)(n+16) ' 

M _ (n+5)(85 n 4 +1817n 3 +11380 n 2 +7072 n-93460) 
1Vl12 ~ 7(n+7)(n+13)(n+15)(n+16)(n+17) ■ 



M 2 1 = - 



(n+6) 2 (170 n 4 +4735 n 3 +42068 n 2 +99767 n-168628) 



12(n+7)(n+12)(n+14)(n+15)(n+16)(n+17) 



and finally 



M _ 4327 n 7 +163698 n 6 +2480127 n 5 +19091004 n 4 +78090428 n 3 +163454544 n 2 
1V±22 14(n+7)(n+12)(n+13)(n+15)(n+16)(n+17)(n+18) 

_l 172290528 re+132098688 



14(n+7)(n+12)(n+13)(n+15)(n+16)(n+17)(n+18) ' 

The matrix 
with 

N _ 2(193 n 5 +5832 n 4 +65284 n 3 +328884 n 2 +727621 n+634422) 



JVll #12 
N 2 1 N 22 



7(n+8)(n+13)(n+14)(n+16)(n+17) 

_ 171 n 5 +47! 

12 8(n+8)(n+14)(n+15)(n+17)(n+18) 
at _ 171n 6 +7C 

iV2i 7(n+8)(n+13)(n+14)(n+16)(n+17)(n+18) 

and finally 



^y i2 _ 171 n 5 +4729 n 4 +45764 n 3 +188570 n 2 +442336 n+1133640 

171 n 6 +7071 ra 5 +H6213 n 4 +959879 n 3 +4245034 ra 2 +10640548 n+15755112 



at _ 4269 n 7 +169934 n 6 +2677678 n 5 +21066480 n 4 +85737209 n 3 +169428298 n 2 
iV22 8(n+8)(n+13)(n+14)(n+15)(n+17)(n+18)(n+19) 

129986220 n-46794888 

8(n+8)(n+13)(n+14)(n+15)(n+17)(n+18)(n+19) ' 

The matrix 

A 7 = 

with 



Pu P12 
P21 P22 



p _ 3(n+5)(129 n 4 +3710 n 3 +36430 n 2 +129960 n+76536) 
ri1 ~ 8(n+9)(n+14)(n+15)(n+16)(n+18) ' 

p _ (n+5)(n+10)(57 n 3 +917 n 2 +2274 ra-11268) 

^ 12 ~~ 3(n+9)(n+15)(n+16)(n+17)(n+19) ' 

p _ -3(57 n 6 +2505 n 5 +44489 n 4 +389955 n 3 +1576582 n 2 +1465908 n-4434696) 

- r2i ~~ 8(n+9)(n+14)(n+15)(n+16)(n+18)(n+19) ' 

and finally 

p _ 2(ra+10)(829 n 6 +27979 n 5 +352571 n 4 +2024521 n 3 +5197384 n 2 +5712396 n+5004720) 
22 ~~ 3(n+9)(n+14)(n+15)(n+16)(n+17)(n+19)(n+20) 

The matrix 

, Qn Q12 

A, 

with 



Q21 Q 



22 



n _ 5(n+5)(n+6)(21n 2 +401n+1920) 
" n ~ 6(n+15)(n+16)(n+17)(n+18) ' 
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Ql2 = 

Q21 = 

Q22 = 

The matrix Ag 



15(n+5)(n+6)(n+8)(n+ll) 
2(n+16)(n+17)(rt+18)(ri+19) ' 

5(n+6)(10 n 4 +329 n 3 +4942 n 2 +36611 ra+96300) 
6(n+15)(n+16)(n+17)(n+18)(n+20) ' 

3(n+6)(ra+ll)(430 n 4 +9773 n 3 +67728 n 2 +129129 n-59220) 
4(n+15) (n+16) (n+17) (n+18) (n+19) (n+21) 







99(n+4)(n+6)(n+7)(n+10) 165(n+4)(n+6)(n+7)(n+8)(ra+10)(n+12) 
K 4(n+16)(n+17)(n+18)(n+19)(n+20) 2(n+16)(n+17)(n+18)(n+19)(n+20)(n+21) ^ 

Notice that if we concentrate our attention on the coefficients within the 
traditional range we see that the first matrix A4 has its first hook made up of 
positive entries, the second hook (which in this example consists of only one 
entry) has negative signs. The second matrix A$ has its first hook negative, 
the second hook positive. The third matrix Aq repeats the behavior of the 
first one, the fourth one Aj imitates the second one, and so on. 

Extensive experimentation shows that this double alternating property 
holds for values of I greater than zero. For coefficient matrices in the tradi- 
tional expansion range, the first matrix has its first hook positive, the second 
one negative, the third one positive, etc. The second matrix has the same 
alternating pattern of signs for the hooks but its first hook is negative. The 
third matrix imitates the first, etc. 

The following picture captures the phenomenon described above for n 
larger than one and when the index k is in the traditional range. 

+ + + ...+ 

+ --...- - + + ...+ 

+ - + ...+ - + -...- 

+ - + - + - 



+ - + - + 



etc. 



5. The relation with matrix valued orthogonal polynomials 

We close the paper remarking, once again, that our matrix valued spher- 
ical functions are orthogonal with respect to a nice inner product and have 
polynomial entries. Yet, they do not fit directly into the existing theory of 
matrix valued orthogonal polynomials as given for instance in Q and 

It is however possible to establish such a connection: define the matrix 
valued function *$>(j, t) by means of the relation 

*(i,t)=$(j J t)$- 1 (0,t). 

It is now a direct consequence of the definitions that the family $>(j, t) 
satisfies all the standard requirements in pi] and not only satisfies a three 
term recursion relation but also *ff(j, i)* satisfies a fixed differential equation 
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with matrix coefficients and only the "eigenvalue matrix" depends on j. In 
other words the family fy(j,t) meets all the conditions given at the begin- 
ning of section 3 and meets also the conditions of the standard theory in |?J 
giving an example of a classical family of matrix valued orthogonal polyno- 
mials. In particular, the coeffients in the differential operator D (obtained 
by conjugation from the one in [jll|] ) are matrix polynomials of degree going 
with the order of differentiation. For a nice introduction to this circle of 
ideas, see the pioneering work in ||. 

In conclusion, we are really indebted to the editors for suggesting a num- 
ber of places where the exposition could be made a bit less terse. One of us, 
A.G., acknowledges a useful conversation with A. Duran that steered him 
in the direction to section 5 above. 
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